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Spin mixing in colliding spinor condensates: formation of an ef- 
fective barrier 



oo 
o 
o 

(N 
> 

o 

OS 



> 
in 



M. Guilleumas 1 , B. Julia-Di'az 1 , J. Mur-Petit 2 , A. Polls 1 



Departament d'Estructura i Constituents de la Materia, Universitat de Barcelona, E-08028 Barcelona, Spain 
Department of Physics and Astronomy, University College London, London WC1E 6BT, United Kingdom. 



PACS 
PACS 



PACS 



03 . 75 . Mn - Multicomponent condensates; spinor condensates 

03.75.Kk - Dynamic properties of condensates; collective and hydrodynamic excitations, 
superfluid flow 

03 . 75 . Lm - Tunneling, Josephson effect, Bose-Einstein condensates in periodic potentials, 
solitons, vortices, and topological excitations 



Abstract. - The dynamics of F = 1 spinor condensates initially prepared in a double- well 
potential is studied in the mean field approach. It is shown that a small seed of m = atoms on 
a system with initially well separated m = 1 and m = — 1 condensates has a dramatic effect on 
their mixing dynamics, acting as an effective barrier for a remarkably long time. We show that 
this effect is due to the spinor character of the system, and provides an observable example of the 
interplay between the internal spin dynamics and the macroscopic evolution of the magnetization 
in a spinor Bose-Einstein condensate. 
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Introduction. — Ultracold atoms trapped by opti- 
cal means are suitable systems to address a broad range 
of problems related to magnetic ordering and dynamics. 
Soon after the pioneering experiments with F = 1 23 Na 
at MIT [1] a number of groups have managed to observe 
spinor dynamics in a variety of conditions, e.g. in a quasi- 
1D system [2] or a 2D trap [3], thus being able to experi- 
mentally address many of the questions posed by theoreti- 
cians [4,5]. 

Among the latter, the existing relation between the in- 
ternal dynamics and the spatial ordering in Bose-Einstein 
Condensates (BEC) provides a beautiful example of how 
the microscopic interactions shape the macroscopic prop- 
erties of the BEC. As discussed recently in [6-8] one such 
example is the formation of magnetic domains in F = 1 
spinor systems (see also [9]). In Ref. [8] the dynamical 
evolution of a ID confined BEC with non-equilibrium ini- 
tial populations in the spin components was studied in 
the mean field approach. As the main emphasis of that 
study was on domain formation, the initial configurations 
considered were always such that there was no relative mo- 
tion between the center of mass of the different Zeeman 
components. 

In this paper we consider the complementary case of 
initial conditions where the spin components are spatially 
well separated inside a harmonic trap, ensuring a cer- 



tain amount of relative center of mass motion between 
them. The subject of collisions between scalar BECs has 
been addressed in several experiments, revealing interest- 
ing features of collective modes [10] and even realizing an 
atom analogue of the Hanbury Brown-Twiss effect of op- 
tics [11]. Other experimental studies on the dynamics of 
spinor BECs [12], as well as binary mixtures [13,14], have 
focused on the evolution of the spatial distribution of the 
different spin components in the system. Specially inter- 
esting is the observation by Hall et al. [13] of a fast con- 
vergence to a configuration with spin-segregation in the 
case of two interacting initially overlapping BECs of 87 Rb 
in different hyperfine states. 

The setup considered here can be seen as a low energy 
confined version of those in Ref. [15], where collisions of 
unconfincd BECs arc analyzed. In the present work the 
system is confined and thus no atoms escape the trap. 

The paper is organized as follows. First the simplest 
extension of [8, 16] is considered, namely the formation of 
magnetic domains in a trapped F = 1 BEC with no rel- 
ative center of mass momenta between the Zeeman com- 
ponents even when the center of mass momentum of the 
system is not zero. Under these premises a full decoupling 
of the Kohn mode [17], dipolar oscillation in the confining 
trap, and the internal spin dynamics is found. 

In the second step we introduce initial spatial separa- 
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tion between the different components, e.g. m = —1 on 
the right side of the trap and m = +1 on the left side, 
which are initially kept separated by means of a gaussian 
barrier. Switching off the barrier at t — the condensates 
collide and the time evolution of the different populations 
is analyzed. 

With these conditions we find a remarkable effect. 
Namely that on F = 1 BECs with initially separated 
to = — 1 and m = +1 components, a small amount of 
to = component produces a barrier-like effect which pre- 
vents the mixing of the to = ±1 components for a long 
time. We further show that this effect is characteristic of 
spinor BECs. in contrast e.g. to binary mixtures, and oc- 
curs in a time scale reachable in current experiments with 
ultracold atoms. 

Description of the system. — In the mean-field 
approach, the F = 1 spinor condensate is described by 
a vector order parameter \1/ whose components ip m cor- 
respond to the wave function of each magnetic sublevel 
\F = l,m) = \m) with to = 1,0,-1. In absence of an 
external magnetic field and at zero temperature the spin 
dynamics of this system confined in an external potential, 
V e:x ^, is described by the following coupled equations for 
the spin components [4,8]: 

i%dil)±i/dt = [H s + ca(n±i + no — n T i)] tp±i 

+C2 V^l , (1) 

ihdipo/dt = [H s + c 2 (ni + n_i)] ipo 

+Pa2^V-i. ( 2 ) 

with H s = -ti 2 /(2M) V 2 + V ext + c n being the spin- 
independent part of the Hamiltonian. The density of the 
m-th component is given by n m (r) = \i/j m (r)\ 2 , while 
n ( r ) = Em IV'mC 1 ")! 2 is the total density normalized to 
the total number of atoms N. The population of each hy- 
perfine state is N m = J dr\ip m (r)\ 2 . Defining the relative 
populations A m = N m /N, it follows that Ai + Ao + A_i = 1 
and the magnetization of the system is M. = X± — A_i. 
The total number of atoms and the magnetization are 
both conserved quantities [8]. The couplings are cq = 
47rft 2 (a + 2a 2 )/(3M) and c 2 = inh 2 (a2-ao)/(3M), where 
M is the atomic mass and ao, & 2 are the scattering lengths 
describing binary elastic collisions in the channels of total 
spin and 2, respectively. The interatomic interactions 
permit the transfer of population between the different 
Zceman components by processes that conserve the total 
spin, |0} + |0) + + 

Preparation of the system. — We consider N = 
20000 atoms of spin-1 87 Rb in a highly elongated trap with 
luj_ = 2-k x 891 Hz and u z = 2tt x 21 Hz [2]. The scattering 
lengths arc ao = 101. 9as and a 2 = 100. 4as corresponding 
to a ferromagnetic behavior, c 2 < [4]. Since uj±_ 3> lo z 
the dynamics takes place along the axial direction and the 
equations of motion become one-dimensional for the lon- 
gitudinal wave functions ip m (z) by rescaling the coupling 
constants cq and c 2 by a factor 1/(2™^), with a± the 



transverse oscillator length [16, 18]. Overlapped with the 
harmonic potential we consider a gaussian barrier which 
separates the system in two symmetric wells (labeled as 
R and L, corresponding to the right and left side respec- 
tively). The resulting confining potential is the same for 
the three Zeeman components, and reads: 

^ext = f^ 2 + ^exp(-z 2 /^ 2 )- (3) 

Given an initial configuration of the system (the relative 
population of spin components inside each potential well: 
X° m , with m = 1,0,-1, j = L.R and A m = A^ + A^), 
the initial state of the system is prepared as follows: (i) 
First the ground-state wave function of the scalar conden- 
sate is calculated from H s on the whole spatial domain. 
The presence of a wide and large enough barrier at the 
center of the trap ensures no overlap between the wave 
function at each side of the barrier. Notice that due to 
the presence of the initial barrier the two wells are not 
parabolic and the wave function in each well is strongly 
non-parabolic, see inset of Fig 2. (ii) Afterwards, the wave 
function of each magnetic sublevel inside each well is ob- 
tained by normalizing the scalar wave function to the de- 
sired initial relative population A^ n in that well [16]. (iii) 
At t = the potential barrier is instantaneously switched 
off and the spinor system is allowed to evolve inside the 
harmonic ID-potential. The dynamics follows by solving 
the system of coupled equations (1,2). Our numerical pro- 
cedure for the time evolution combines the split operator 
method with the fast Fourier transform to treat the ki- 
netic terms and a fourth-order Runge-Kutta method for 
the remaining terms of the dynamical equations [8] . 

Scalar condensate. — In order to understand the 
spinor-driven effects on the condensate dynamics, let us 
first address a spin-polarized condensate of to = —1 atoms 
all of them localized on the right side of the trap: A^=0 
(for all to), = Xq = and A^ = 1. Since all the atoms 
are in the same hyperfine state this case is equivalent to 
having a scalar condensate. At t = the gaussian barrier 
is switched off and the system evolves inside the harmonic 
trap bouncing from right to left in the well-known Kohn 
mode [17]. 

In Fig. 1 we depict the fraction of atoms in the left side 
of the trap as a function of time (solid line) . For compari- 
son, the dipole oscillation of the ground state of a displaced 
harmonic trap is also displayed (dashed line). The Kohn 
mode is clearly identified by the frequency of the popula- 
tion oscillations being equal to the axial frequency of the 
trap ui z . Further, as our initial configuration is not the 
ground state of the harmonic trap, the Kohn oscillation 
is convoluted with other frequencies associated with the 
deformation of the density profile that bounces inside the 
trap. 

Component driven dynamics. We proceed now 
to introduce spin dynamics in the system. To this end, 
we analyze the time evolution of a system with all three 
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Fig. 1: Fraction of atoms on the left side of the trap as a 
function of time for two different initial states of a scalar BEC 
(sketched in the inset): (a) harmonic trap with a gaussian bar- 
rier (solid line) and (b) ground state of a displaced harmonic 
trap (dashed). 

Zeeman components populated, and the atoms initially in 
the right side of the trap. 

If we start with Af x = 25% and A^ = 50% (and 
A^ = Vm), which corresponds to the ground-state con- 
figuration in spin space [8] , we observe no transfer of pop- 
ulation among the different components, while the total 
density follows the same behaviour as the solid line in 
Fig. 1 (Kohn mode). Moreover, if we start from a spin 
configuration away from the ground state, the total den- 
sity still follows a dipole-like oscillation, on top of which 
a transfer of population happens identical to that occur- 
ring in the case of a system at rest in a simple harmonic 
trap [8, 19]. In other words, the center-of-mass motion 
and the internal (spin) dynamics are fully decoupled, as 
expected for the dipole mode in a harmonic trap [20] . We 
remark that the time evolution of X m (t) is exactly the 
same as in the system at rest. The relative populations 
oscillate around the corresponding equilibrium values at 
a given magnetization (e.g., Ao = 50%, A±i = 25% for 
M. = 0), as was already pointed out in Ref. [16]. 

We have seen how the presence of internal spin dynamics 
does not affect the evolution of the total density n(z, t) 
of the system when all the atoms have the same initial 
spatial distribution. Now we will illustrate how the exact 
initial conditions of the system (spin configuration and 
spatial distribution thereof) can have a strong effect on 
the internal and the spatial dynamics, showing genuine 
effects associated to the spinor nature of the condensate. 

Let us consider the case where initially all m = 1 
(m = —1) atoms are located on the left (right) side of 
the trap, with a vanishing overlap between the two con- 
densates (f_ g dz[n_i(z) + n + i(z)\ = for 5 = 4/im). A 
small amount of atoms with m = is added on both sides, 
A o = = 1%, together with Xf = A^ = 49% and 
A_i = Af = (see inset of Fig. 2). This small admixture 
of m = allows the exchange of atoms between spin com- 
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Fig. 2: Time dependence of the total (top) and fraction on the 
left side (bottom) relative populations of the various Zeeman 
components of a spinor condensate, for the initial conditions 
Ai = \ R 1 = 49%, \o = \q = 1%. Solid, dotted and dashed 
correspond torn = — 1, m = and m = +1 respectively. Notice 
that the lines for m — ±1 in the top panel are indistinguishable. 

ponents: |0) + |0) <-> | + 1) + | — 1) within the approach of 
Eqs. (1,2). Then, the population of each sublevel changes 
with time but due to the symmetry of the equations and 
the initial conditions we have that for any t, X\ (t) = A_i (i) 
(see Fig. 2). 

The resulting dynamics after switching off the trap is 
as follows. At short times, m — atoms are created at 
the center of the trap where the overlap between m = ±1 
components is larger [cf. Eq. (2)]. Here the m = ±1 com- 
ponents fuse together giving rise to m = atoms with a 
conversion rate proportional to c<x. This can be seen in 
the lower panel of Fig. 2, where we present the relative 
populations on the left side of the trap as a function of 
time. We observe after ~ 5 ms a temporary saturation of 
the number of m = — 1 atoms on the left side for up to 
350 ms, while the m = +1 component keeps losing pop- 
ulation, which is transferred to the m = Zeeman state. 
The oscillations seen in all X J m in Fig. 2 are reminiscent 
of the Kohn mode discussed above, and have a dominant 
frequency close to that of the harmonic trap (u z ). 

The long time behaviour observed in Fig. 2 is deter- 
mined by the initial conditions that correspond to a highly 
excited configuration. After the clouds move towards each 
other and collide, the initial collision energy is distributed 
among internal degrees of freedom, thus acting as an ef- 
fective temperature. Therefore the populations perform 
damped oscillations around the equipartition configura- 
tion [16]. 



M. Guillcumas, B. Julia-Diaz, J. Mur-Petit, A. Polls 




m=+l 



-25 25 50 
z (|Lim) 



Fig. 3: Density distribution of atoms in the trap, at time t = 
115 ms with Ao(i = 0) = 2% (top) and for the case of Xo(t = 
0) = (bottom). Solid, empty and dashed correspond to m = 
— 1, m = and m = +1. A ~ 4/im spatial binning has been 
employed to generate the plot. 



To further understand the effect of the to = compo- 
nent, we show in the upper panel of Fig. 3 the density 
profile of the various Zeeman components at t = 115 ms, 
well into the time evolution of the system. We observe 
the presence of to = atoms located at the center of the 
trap, with m = 1 (to = —1) atoms remaining to its left 
(right). In summary, the initial to = ±1 components are 
steadily converted into m = atoms which stay mostly 
at the center of the trap, and act as an effective barrier 
that prevents the spatial mixing of the surviving m = ±1 
components. These stay mostly on their original side of 
the trap for several hundreds of miliseconds, with only a 
small fraction of atoms crossing over to the other side. 
The long lifetime that the to = barrier shows in this 
case is to be contrasted with the much faster spin dynamics 
seen previously, where the typical timescale for population 
transfer between initially overlapping Zeeman states was 
~ 100 ms [8]. It also contrasts with the spin-segregation 
time scale observed in binary mixtures of 87 Rb [13], as will 
be addressed below. Interestingly, qualitatively similar re- 
sults have been obtained in an antifcrromagnctic system 
simulated by changing the sign of C2 (pz > 0). 

We have also observed that even a very small m = 
seed in one side (0.02%) has the same effect of acting as 
a barrier and stops the mixing of the m = ±1 atoms. 
However, as expected, different Xo(t = 0) result in different 
transient behaviours. In particular, the larger the to = 
seed the longer the time necessary to have a mixture of 
the other Zeeman components with similar concentrations 
in both sides. Nevertheless, the long time (~ 400 ms) 
configuration is roughly the same. 

To emphasize the spinorial origin of this barrier effect we 
compare with the mixing process of two interacting BECs. 
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Fig. 4: Fraction of atoms in the m = — 1 (solid) and m = +1 
(dashed) Zeeman component on the left side of the trap, as a 
function of time for a binary mixture. The inset shows the ini- 
tial spatial distribution together with the confining potential 
before switching off the a = 5/xm gaussian barrier. The dot- 
ted line stands for the time evolution of A_i(t) for an initial 
configuration given by a = 40/^m. 



Such binary mixtures have been realized experimentally 
by coupling two Zeeman states via a two-photon drive. 
For example, mixtures of 87 Rb in the hyperfine states \F = 
2,m = 2) and |2, 1) [10] as well as |1,-1) and |2, 1) [13,14] 
have been both produced. In the mean-field framework 
their dynamics can be easily simulated with the coupled 
differential equations (1,2) by taking a zero amount of m = 
atoms. 

The preparation of the system follows the same steps as 
above, setting now the initial conditions as = A_ x = 
50% and A+ x = A^ x =0. The barrier is switched off at 
t = and the two condensates move inside the harmonic 
trap towards each other, colliding after ~ 10 ms. As can 
be seen in Fig. 4, at this time a sizeable amount of atoms 
of each component has already reached the other side of 
the trap, and rapidly mixes with the other component. 
As expected, A_j grows initially with a similar rate as it 
did in Fig. 2, but departs from that behaviour after times 
~ 5 ms, due to atoms coming from the other side of the 
trap. There is also a similar decrease of X I f_ 1 since an 
equal amount of to = 1 atoms goes to the right. Around 
t = 25 ms the two components are mixed and there is 
almost the same number of to = ±1 atoms in each side of 
the trap, i.e. there is a fast and complete mixing of both 
components, as shown in the lower panel of Fig. 3. 

If the width of the initial gaussian barrier is increased, 
the initial collision energy is larger, and thus the number of 
m = — 1 atoms reaching the left side is also larger initially. 
However, it reaches the same asymptotic value as before 
(cf. dotted line in Fig. 4). Therefore, we conclude that the 
transition to a steady state is much slower in the spinor 
case than it is for a binary mixture (cf. Figs. 2 and 4), 
which qualitatively agrees with the observations in [12,13]. 

Conclusions. — In summary, we have studied the in- 
terplay between the internal spin dynamics and the spatial 
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evolution of a spin-1 condensate in a confining potential. 
We have described the effect produced by the presence of 
m = atoms in a system initially prepared with two spa- 
tially separated spin components m = +1 and m = — 1 . 
An initial small population of the m = component dras- 
tically affects the mixing time of the system as compared 
to a binary mixture, producing a barrier-like behaviour 
which keeps the other two components from mixing for 
times up to ~ 350 ms which is a timescale on which ex- 
perimental observations should be feasible. We have also 
demonstrated the genuine spinor character of this effect, 
and that it does not appear in scalar condensates or bi- 
nary mixtures. The population transfer term of Eqs (f ,2) 
plays a key role in the generation of the effective barrier. 

The seed of m = atoms can be a residual of the prepa- 
ration of the initial state in an experiment. A similar effect 
would also be produced in an initially pure, | + f) |— 1) mix- 
ture by two-body spin- flip processes | + 1) + | — 1) — » |0) + |0) 
at the initial stages of the evolution. Therefore, the pres- 
ence of the effective barrier preventing the mixing of spin 
components appears as a general phenomenon in spinor 
BECs prepared with spatially separated Zeeman compo- 
nents. 

fn the future, we will address more general examples of 
interplay of spin and spatial degrees of freedom in spinor 
Joscphson junctions whose barriers have a more complex 
time dependence, and can distinguish between different 
hyperfine components [21]. We note also that more ex- 
otic spatial orderings and dynamics are expected to occur 
in optical lattices, where oscillation in the populations of 
two different Zeeman components have already been re- 
ported [22]. 
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